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INTRODUCTION 

This study describes a method for obtaining an 
approximation to the auto- and cross-correlation 
functions, R X x anc * R xy. Tne method involves use 
of the EAI PACE® TR-48 General Purpose Analog 
Computer with Electronic Mode Control, used with 
the EAI DES-30 Digital Expansion System ... a 
low-cost, general-purpose, patchable digital logic 
system. 

This logic package, although capable of operating 
autonomously as an aid in digital design or instruc- 
tion, was designed primarily as an expansion to 
the TR-48 computer to provide basic hybrid 
capabilities for the small computer facility. Thus, 
the necessary logic and control functions for the 
hybrid implementation of this application are pro- 
vided by the digital capabilities of this compact, 
completely self-contained unit. 

Specifically, the study describes the implementa- 
tion of an approximation to the auto-correlation 
function 

T 

R^ = lim Jjr J x(t) x(t - r) dt 

The extension to the cross-correlation function, 
R xv , is obvious and requires no changes in the 
program. 

METHOD OF SOLUTION 

The problem is to plot R^ (r). A train of blips 
of period f (variable) is generated. On each such 
blip the product x(t) x(t-r) is sampled and then 
held on one track/ store unit; then, the present 
value of x(t) is also sampled and held... it is used 
to represent x (t-r) on the next sampling pulse. 
Now the entire affair is run in the RO mode, with 
t changing very slowly (so it may be considered 
constant for at least one RO cycle). A standard 
tandem track/ store arrangement picks up and 



stores the present value of Rxx for presentation 
to the plotter. 

THE PROGRAM 

The program is shown in Figure 1. The compo- 
nents are labeled sequentially and do not represent 
an actual patching. 

DA switch amplifier 08, multiplier 10, integrator 
11 and electronic comparator 12 comprise the 
variable frequency pulse generator which gener- 
ates blips with period T . This is a standard 
HYDAC program. The blips are obtained at the 
output of differentiator 07. 

MT06 generates two consecutive 100 /lis signals. 
It does so as follows; the timing blip sets the MT 
through OR gate 17A. It also sets FF26A, which 
enables AND gate 17B. When the MT returns to 0, 
DIF 06 produces a blip which sets MT06 again for 
another 100 jLts period (DES-30 monos are 100% 
duty cycle). The same blip from differentiator 06 
also resets FF26A so it is ready for the next 
timing blip. Note that this means that r cannot be 
less than 200 jus, and a practical minimum value 
is probably 400 as. The timing just described is 
shown in Figure 2. 

The output of MT 06 is 1 'time- shared": When 
FF26A is set, the first 100 j^s signal is routed 
through AND gate 17E to T/S 04 which tracks and 
holds the product of the present value of x, x(t), 
and the 1 'delayed' ' value x(t-r). The second 100 j^s 
signal occurs after FF26A is reset , and it is 
routed through AND gate 17F to T/S 00 which 
tracks and holds x(t) (for use as x(t-r) on the 
next timing blip) . 

Integrator 05 performs the integration involved in 
generating R„ The pot preceding integrator 05 is 
set to Ki/2T where is a scale factor and T 
(not to be confused with r ) is the RO compute 
period, or the limits of the integration. 
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+ X(t) X(t-T ) 




IC (RS) 



Figure 1: Computer Program 



CLOCK 

TIMING SLIP 
MT 06 

AG 100" 
AG 101 



JJJJ LL 



1111 LL 



MM LL 



H — 

TL 



J~L 



U 1 I U L 



N- 100-H N-IOO-rt 
/xS fiS 



J L 



J L 



r — l 



J LT~ 

J br 



Figure 2: Timing Diagram ( r constant) 



Track/ stores 06, 07 are a standard tandem scheme 
for plotting a "high speed' ' variable in slow speed. 
The value of at the end of each RO cycle is 
determined by T/S 05 and transferred to T/S 06 
where it is held unchanged until updated after the 
next run. The signals to operate T/S 06, 07 can be 
tapped from any integrator mode control bottle 
plug. 

Integrator 13 controls the parameter r and also 
provides a drive proportional to r to the plotter 
arm. To obtain a "normal" plot (increasing r), 
integrator 13 should start from maximum (+ 10V) 



and ' 'integrate down" to 0. This corresponds to 
starting from the maximum pulse rate at DIF 07 
(minimum r) and decreasing frequency (increasing 
T). The plotter arm is driven from -10 volts to 0, 
and proper parallax adjustment will have to be made. 

The minimum value of r can be varied by changing 
the feedback capacitor and/or changing the gain of 
integrator 11. Note that all capacitor values are 
obtainable by playing with the 0.1/3 and NORMAL/ 
SPECIAL bottle plug. Table 1 summarizes the 
situation: 



TABLE 1, MINIMUM r PERIOD 



Normal 


Special 


0.1/5 


Capacitor 


Gain 1 


Gain 10 


OUT 


OUT 


IN 


10 fxY 


2 sec 


200 ms 


OUT 


OUT 


OUT 




200 ms 


20 ms 


OUT 


IN 


IN 


.02 /xF 


4 ms 


400 


OUT 


IN 


OUT 


. 002 M F 


400 M s 


40 ms 



As indicated previously, 200 ijls is the lowest 
practical r since the track/ stores 00"," 03 require 
each 100 /xs to track properly. 

Note that integrator 11 is always in OP; integrator 
13 operates at a normal rate, while integrator 05 
is in the RO mode. Also note that the output of 
integrator 13 must not be allowed to go negative: 
the operator should monitor its output and flip 
FS 1 (also lift pen!) to the open position when the 
output gets close to zero. An automatic program 
to do this can easily be implemented but, unless 
plotter pen control is also available, the value of 
this addition is not justified. 

LIMITATIONS 

The substitution of track/store units for a con- 
tinuous delay device imposes a severe limitation 
on T. The value of r may not exceed 1/100 of a 
cycle of x(t), or about .07 radians. Placing a 
linear interpolator between track/ store 04 and in- 
tegrator 05 may relax this limit somewhat. At .07 
radians, the error is probably of the order of 1%. 



A NOTE ON SYMBOLS 

In this memo the following symbols are used: 



TRACK/STORE 



DA SWITCHES 



ELECTRONIC 
COMPARATOR 




digY 



BUFFERED (CLOCKED) 
A->D TRUNK 



Figure 3: Definition of Symbols 

All other symbols used are standard analog or 
HYDAC symbols. 
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